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SUMMARY 

A counter rotating bearing rig has been designed and constructed to study 
transient elastohydrodynamic lubrication phenomena. This paper describes new 
instrumentation and documents test procedures. Ball and race speed measure- 
ment systems and the capacitance (film thickness) measurement system were 
upgraded. Methods for measuring bearing torque and race temperatures were 
implemented. 


INTRODUCTION 

Starvation theory, a subdivision of elastohydrodynamic lubrication (EHL), 
was described a number of years ago (ref. 1). It relates to the situation 
occurring in ball bearings having a restricted oil supply in which pressure 
build-up in the inlet is inhibited resulting in a film thickness thinner than 
calculated by classical EHL theory (ref, 2). 

Starvation theory fails to adequately describe instrument ball bearing 
behavior. Another subdivision of EHL, "parched" describes this behavior where 
there is no free bulk oil in the system (ref. 3). The lubricant films are so 
thin that they are immobile outside the Hertzian contact zone. This regime is 
of practical importance since parched bearings require the least driving 
torque and also have the most precisely defined spin axis. 

The lubricant inside a Hertzian contact is not immobile and is eventually 
squeezed out resulting in long-term film thinning transients (refs. 4 and 5), 
which leads to film failure and subsequently bearing failure. In order to 
prevent failure, make-up oil must be added to balance the loss due to cross 
flow. 


As film thinning occurs, shear rate increases and mechanical shear energy 
is then concentrated in smaller and smaller Hertzian volumes. At some point 
in this process, the shear energy causes irreversible chemical changes in the 
lubricant, resulting in degradation and/or polymerization (ref. 6). 

For retainerless bearings operating at low contact stresses, it has been 
suggested that a beneficial in-situ polymeric film is produced which yields a 



geometrically perfect, plastic ball separator in just the right location 
(ref. 7). This is particularly important since elimination of the retainer in 
a bearing precludes any retainer instability problems which commonly occur in 
gyroscope bearings (ref. 8). 

To study these concepts, a Transient Elastohydrodynamic Lubrication 
Apparatus (TEHL) was initially designed and constructed by Kingsbury (ref. 3). 
Later, modifications by Hunter (ref. 9) allowed for calculation of film thick- 
ness from capacitance measurements. 

The objective of this paper is to describe several new modifications to 
the TEHL apparatus and to completely document its operation. The ball and 
race speed measurement systems and the capacitance measurement system were 
upgraded. Methods for measuring bearing torque and race temperature were 
implemented. Finally, a computer data acquisition system was installed. 


TEHL APPARATUS - BASIC OPERATION 


The overall apparatus is shown in figure 1. It consists of the TEHL 
apparatus itself, inner and outer race drive motors (the latter containing a 
torque sensor) and a loading mechanism. The TEHL apparatus appears in 
figure 2. The upper bearing is the test bearing. The outer race of this 
bearing is driven by a synchronous hysteresis motor through a toothed belt 
drive. A keyed bushing is press fit into the inner race. The load ring is 
keyed into the bushing and a spindle, which is driven by a second synchronous 
motor through another toothed belt drive. This allows the bearing to be run 
in a counter - rotating mode (the races spin in opposite directions such that 
the ball complement is stationary). The bearing is loaded with a deadweight 
axially through the load ring and the central load shaft that extends down 
through the spindle (fig. 3). The four lower bearings provide support and 
alignment . 


Ball-race coupling is followed by measuring the basic speed ratio ( BSR) 
of the bearing (ref. 10). BSR is defined as: 


BSR 



( 1 ) 


where Q is the ball spin rate and H and fl are the inner and outer 

b *■ i o 

race spin rates, respectively (algebraic sign defined as shown in fig. 4). 


BEARING CAPACITANCE, FILM THICKNESS, AND CONDUCTIVITY 

The bearing capacitance measurement system consists of a Hewlett-Packard 
4267A (Inductance-capacitance-impedance) LCZ meter, interfaced with an IBM 
compatible personal computer through an IEEE-488 interface bus. The LCZ meter 
is connected to the bearing by four coaxial cables. Two of these cables sup- 
ply a known alternating current across the bearing and two carry the voltage 
that is developed across the bearing as a result of this current back to the 
LCZ meter as shown in figure 5. Measurements of the amplitude and phase shift 
(relative to the supply current) of the returned voltage are made by the 
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instrument , and the total impedance across the bearing is calculated* The 
impedance is resolved into an effective parallel capacitance and conductance, 
which are the actual values displayed by the meter. The average contact film 
thickness is then calculated from the capacitance readings by a computer 
program from a method developed by Dyson et al. (ref. 12) for cylinders but 
successfully applied to ball bearings by Allen, et al. (ref. 13). 

Actual electrical connections to the bearings are made through two sets 
of slip rings (fig. 5). The first set carries all four wires from the 
stationary lab reference frame to the rotating outer housing. Two wires (one 
current supply and one voltage sense) are attached to a screw, which in turn, 
is tightened against the outer bearing race. The other two wires pass through 
the second slip ring assembly (whose stator rotates with the outer housing) to 
the load ring, where they are attached by a screw. The load ring makes direct 
metal contact with the inner race, completing the circuit. The whole appara- 
tus is electrically isolated from ground (except for a 5 Mfl bleed-off resis- 
tor, which drains away static build-up). Isolation from ground is necessary 
for proper operation of the LCZ meter. Typical capacitance readings are on 
the order of 500 pF, with an accuracy of about ±2 pF, and conductance 
readings are in the 1 ji s (microsiemens) range, accurate to ±0.1 /is. 

Stray capacitance and conductance values are measured by breaking the 
circuit and by backing out the outer race contact screw, and then taking 
readings from the meter. These stray background readings are then subtracted 
from the values read from the LCZ meter. 


BALL AND RACE SPEEDS 

The ball and race speed measuring circuit is more complicated than the 
capacitance measuring circuit (fig. 6). It is also considerably more accu- 
rate, on the order of ±0.001 percent. This is important in the calculation of 
the rate of slip in the ball-race contacts. There are three basic sections to 
this system: the race speed measurement, the ball spin measurement, and the 

ball complement position controller. 


RACE SPEED 

The races are driven by synchronous motors through timing belt and pulley 
arrangements. To determine the race speed, it is sufficient to measure the 
drive frequency of the motors, and divide by a constant (in this case, the 
integer 4). The constant is determined by the pole arrangement of the motors, 
and the ratio of the number of teeth on the driving pulley to the number on 
the driven pulley. The signals from monitor outputs on each of the variable 
frequency power supplies used to drive the motors are fed into a pair of 
counters. The counters measure the period of the signals and transmit the 
information to the computer through the IEEE-488 interface bus. 


BALL SPIN RATE 

To measure the ball spin rate, one of the bearing balls is slightly 
magnetized. A pick-up coil, made of 1000 turns of number 34 A.W.G. magnet 
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wire wound on a laminated iron core, is positioned near the magnetized ball, 
(The ball is stationary in space because of the counter-rotating mode of 
bearing operation, as noted previously,) The output signal from the coil is 
on the order of several millivolts, and has significant components of many 
frequencies related to the apparatus and the instrumentation. This signal is 
first boosted by a preamp with a gain of around 100, and a low pass filter 
with a roll-off to 12 dB/octave past the -3 dB point of 1000 Hz. (Ball spin 
rates are typically 200 to 300 Hz.) The signal is then sent through a series 
of three sharp (Q = 50 to 100) Deliyannis band pass filters (ref. 11), an 
automatic gain control, a differentiator, and then a zero-crossing detector. 
The output is a square wave of precisely the same frequency as the ball spin 
rate. This is then sent to a counter, which measures the period of the wave, 
and sends the information to the computer through the IEEE-488 interface bus, 
just like the race speed counters. 


BALL COMPLEMENT POSITION CONTROLLER 

An important part of the speed measuring circuitry is the ball complement 
position controller. In order for the ball spin rate to be accurately 
measured under the assumption of zero ball complement orbital rate, the 
complement velocity must be exactly zero. The controller makes minute 
adjustments to the inner race speed to keep the ball complement stationary 
during the course of the test run. 

A second coil, identical to the one used in the ball spin rate measure- 
ment circuit, is placed near the bearing, about 45° away from the first coil 
During a test run, the magnetized ball sits between two coils, generating an 
EMF in each of them. If the ball is exactly halfway between them, the 
magnitude of the EMF generated at the ball spin frequency will be the same in 
each coil. If the ball begins to move toward one coil, the EMF generated in 
that coil will increase in magnitude, while that in the other coil will 
decrease. The difference in magnitude gives a measurement of the magnetized 
ball's position between the two coils. 

The second coil is fed through a preamp identical to the one used with 
the first coil. The output signals from the two preamps are then fed into the 
controller box (fig. 6). Inside the box, both signals pass through identical 
parallel Deliyannis filter sections, and then rectifier sections, to obtain 
the two amplitudes. The difference in these two amplitudes is sent through a 
buffer/damper section and then into a PID controller. The output from the 
controller may go through an invertor, if necessary, to get the proper 
controller action, depending on the direction of motor rotation. The output 
is then added to a dc bias voltage to form the main control signal, which is 
output from the controller box. This control signal is sent into the voltage 
controlled oscillator input of a signal generator, which in turn is used to 
synchronize the inner race drive motor power supply. Thus, any disturbances 
in the position of the ball complement results in an adjustment of the inner 
race rotation rate to return the ball complement to its original position. 
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TORQUE AND TEMPERATURE MEASUREMENTS 


In setting up bearing torque measurement system consideration had to be 
given to not disturb the other measuring apparatus. In order to accomplish 
this, it was necessary to use a somewhat unconventional method of torque 
measurement . 

It is relatively easy to measure the torque required by the test bearing 
and the upper support bearing together: Simply measure the torque output of 

the motor, multiply by the ratio of the number of teeth on the driven pulley 
to the number on the driver pulley, and subtract losses due to windage, and 
slip rings. Because of the mechanical coupling between those two bearings, it 
is necessary to consider thermal effects to determine what fraction of the 
total torque is actually required by each bearing. Since both bearings rotate 
at precisely the same speed, the energy dissipated in each bearing is propor- 
tional to the torque generated by each bearing. The constants of proportion- 
ality are the same. 

Figure 7 shows the simplified heat flow diagram of the outer housing. 

The top and bottom of the housing can be treated as insulated surfaces, 
because the plexiglas covers have a much higher thermal resistance than does 
the metal, or thin bakelite insulator. The bearings are also separated by a 
bakelite ring, which prevents large quantities of heat from flowing from one 
bearing directly to the other. As a first order approximation, the housing 
and the two races are assumed to be isothermal because of their extremely low 
thermal resistance. If the housing is thermally symmetric, then the ratio of 
the temperature difference between the test bearing and the housing, and the 
difference between the support bearing and the housing is exactly equal to the 
ratio of the test bearing torque to the support bearing torque. Knowing this 
ratio, and the sum of the torques, the test bearing torque can be determined. 

A more general treatment of this concept is contained in appendix A. 

Motor torque is measured using a reaction torque sensor (fig. 1). The 
motor is suspended from the sensor, so that any torque applied to the drive 
pulley is counteracted by the sensor. The sensor is a strain gage type, so a 
conditioner/indicator is also used. The conditioner/indicator has an analog 
output, which is used to transmit data to the computer. Typical torque values 
are on the order of 0.071 N-m (10 in.-oz) and are accurate to ±0.0007 N-m 
(±0.1 in.-oz) . 

Temperatures are measured using thermistors, since they are more accurate 
than thermocouples at temperatures near room temperature. The thermistor 
resistances are measured using a Fluke 8520A Digital Multimeter ( DMM ) four- 
wire system, which, in principle, operates much like the LCZ meter four-wire 
system (two wires carry a known current, two wires carry back the resulting 
voltage drop) , except that a direct current source is used, instead of the ac 
source used by LCZ meter. The DMM also communicates with the computer using 
the IEEE-488 interface bus. 

In order for a single DMM to read the resistance of all four thermistors 
(test bearing, support bearing, outer housing, and ambient air temperatures), 
the leads have to be switched from one thermistor to the next. This is done 
with a relay box, using gold plated contact relays. It is important to use 
gold plated contact relays in good condition to keep the lead resistance to a 


5 



minimum. Other types of relays, or even gold plated contacts with a little 
wear, produce variations in lead resistance, which show up as noise in the 
temperature data. 

The DMM leads pass through the same slip ring system as the LCZ meter to 
get to the outer housing. Screw terminal connections are used to attach the 
leads from the slip rings to the leads from the thermistors. Because of the 
limited number of terminals on the slip ring assembly, it was necessary to 
arrange the thermistors in a "Y" arrangement. 

The thermal inertia of the bearing housing causes a time lag in the 
response of the temperature readings to variations in torque. It was assumed 
that the temperatures (K ) are related to the bearing torques (T i ) through the 
Laplace domain transformation: 


K 4 (s) ! 

G(s) - = 

T 1 (s) T s + 1 


( 2 ) 


This is known as a first order lag. When a step change in torque occurs for 
some reason, the temperature asymptotically approaches its new value according 
to the governing equation, T * C e" {t/T> , where T is a response time con- 
stant and C is a general constant. It is possible to remove the time lag 
from discrete data by applying the transformation: 




1 - e 


-(t/T) 


( 3 ) 


(See appendix A). Here the Y is the temperature with the time lag removed, 
and t = t - t , where t represents the elapsed time. The subscript i 
represents the value at the ith discrete value, and T is the response time 
constant. This transformation is differential in nature, and thus, any noise 
in the data is magnified. Actual values for T for the various thermistors 
are determined from plots of temperature from runs where the bearing load was 
abruptly changed, creating a step in the bearing torques. 


EXAMPLES OF DATA 

Figures 8 to 10 show examples of data generated with a super-refined gyro 
bearing lubricant (SRG) using the above described instrumentation. Table I 
contains typical physical property data for SRG 200. Table II lists the 
bearing geometry and nominal operating parameters. 

The figures show typical ranges, variations, and resolutions for the 
different measured (conductance, capacitance, race and ball speeds, tempera- 
tures and motor torque) and calculated values (film thickness, basic speed 
ratio, and bearing torques). All plots were generated from a single test that 
was run to failure. Failure occurred at approximately 18 000 sec. 
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CONCLUSIONS 


A fully instrumented and automated bearing test rig now exists that can 
be used in the study of the failure of angular contact ball bearings running 
in the parched EHL regime. This real-life bearing test should give quantita- 
tive results indicating which factors are more significant in the failure 
sequence. This, in turn, should help in the synthesis of new lubricants and 
bearing materials and coatings, and provide valuable data for the design of 
systems which operate in the parched regime. 


APPENDIX A 


The test bearing torque can be viewed as the product of two quantities: 
the sum of the test and support bearing torques, T,, and the fraction of that 


sum which goes to the test bearing, 
expressed as follows: 


13 . Mathematically these quantities are 




T Z - r s + T t 

(1A) 



P* - T t /T* 

(2A) 

where : 



r t 

test bearing torque 



T 

Q 

support bearing torque 



Tjj is found by measuring the 
the ratio of pulley diameters 
subtracting out losses due to 

torque output of the drive motor, multiplying by 
to obtain the torque applied to the housing, and 
slip ring friction and windage: 



D h r m 

r s = - r, 

S D m 

m 

(3A) 


where: 

torque output by the motor 
T l torque lost to windage and slip ring friction 
D h diameter of housing pulley 

D diameter of motor pulley 

Losses due to windage can be shown to be on the order of the resolution of 
the motor torque sensor and are neglected in practice. Loss due to the slip 
ring friction is a constant, experimentally measured to be 0.0044 N-m 
(0.62 in. -oz) . 

To determine 13^, it is useful to define a new quantity 7}^, the ratio of 
the support bearing to the test bearing torque. The ratio is related to B t 
by the expression: 


P t - 


i 

1 + v t 


(4A) 
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Since both bearings rotate at identical speeds, 


_ T s _ _ Q> 

Vt ’ ~t " ’ sk 

where: 

Q rate of heat dissipation by support bearing 
s 

Q t rate of heat dissipation by test bearing 
(] support bearing speed 

s 

test bearing speed 

From this relationship it follows that: 

Qt 
Qt + 


(5A) 


(6A) 


In the determination of ft ’ , it is convenient to think of the housing as a 
network of thermal resistances (fig. Al), and to draw on an analogy between 
heat transfer theory and electrical circuit theory. In this model, the 
thermal resistance of the bearing races and the metallic part of the housing 
are neglected, since their thermal conductances are two orders of magnitude 
higher than the plexiglass or bakelite parts, while the geometric dimensions 
of all the parts are similar. and R 2 represent the thermal resistance 

between the metallic part of the housing and the test and support bearing 

races, respectively. R represents the thermal resistance between the 

bearing races, while R represents the thermal resistance between the 

metallic part of the housing and ambient air. R 3 and R 4 are included in 

the model to account for other minor means of heat dissipation from the 
bearing races. 


Ideally, R 3 = R 2 , and 
and Q are given by 


'4' 


and 


would be infinite. In this case, 


T - T. 


[7A> 


Qt 


T - T 
x t 1 h 


[8A) 


where: 

T temperature of the 
T g temperature of the 
T h temperature of the 


test bearing race 
support bearing race 
metallic part of the housing 
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It follows that is given by 

p t -e ( 9 a ) 

where 8 is the temperature difference ratio defined by: 

T - T 

e - 1 h (ioa) 

T, + T - 2T h 
t 8 n 

Since the idealized case is not realized in practice, the more general 
case of figure Al, with all the R's finite and different is used to solve for 
P t as a function of 8 . A Taylor series expansion of the resulting expres- 
sion yields the following polynomial in 0. 


A 


= A + [1+B+C]‘ 


oo 


£ d n D 


n_n-l 


n=l 


(11A) 


where : 

(X 2 - Xj) ( 1 +8^ 

A * 7 1 

[1 + X 2 ( 2 + e 2 ) j ( 1 + e 3 ) + X 3 (l + e 2 ) 

2Xj(l + e 3 ) - X 2 (2 + 2« 1 + e 2 + e i e 2 ) - e 2 X 3 - 
B * . 

[1 + ^ 2 (2 + e 2 ) J { 1 + e 3 ) + X 3 (l + e 2 ) 

(1 + e 1 )(X 2 - Xj) fejl + X 2 (2 + e 2 ) J ♦ e 2 X 3 } 

C * 1 l 

([1 + X 2 (2 + e 2 ) J ( 1 + e 3 ) + X 3 (l + e 2 )j 2 

ei[i + X 2 (2 + « 2 )J + e 2 X 3 
jl + X 2 ( 2 + e 2 ) (l + e 3 ) + X 3 ( 1 + e 2 ) 


where : 


X^ = R 2 /R x ; X 2 = ^ c /^3 f and X^ ** R^/R 3 

are the dimensionless "leakage parameters,** and: 

e 1 = (R/Rj) - 1 ? and e 2 - (R 3 /R 4 ) - 1 

are the dimensionless "asymmetry parameters." Under the idealistic assump- 
tions given above (R i = R 2 , and R 3 , R 4 , and R x infinite), all the leakage 
parameters and asymmetry parameters are zero, and equation (11A) reduces to 
equation (9A). Under the more realistic set of assumptions that the leakage 
parameters and asymmetry parameters are nonzero, but small compared to unity, 
the following expression is obtained from equation (11A): 
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Pt = < X 2 " X l> 


+ 


(12A) 


1 - 2(X 2 - X x ) - e x p + ifl 2 


As expected, this expression is very similar to equation (9A), except for a 
small correction to the linear term, a small constant, and a small nonlinear 
term, A combination of theoretical and experimental considerations have shown 
that s 

x 2 - X x < 0.03 and |e 1 | < 0.02 

Using these limits, the maximum error in using the simplified equa- 

tion (9A) is less than 9 percent. Therefore, equation (9A) was used for data 
analysis. Though this method is quite effective for steady state or slowly 
changing conditions, the thermal inertia of the housing makes the dynamic 
response of the system rather poor, with time constants on the order of 
10 min. To deal with this problem, three new variables are introduced: , 

Y g , and Y h . These represent the values that the temperature T t , T g , and T h , 
respectively would have if the thermal inertia of the system was negligible. 
Transformed into the Laplace domain, the actual temperatures are assumed to be 
related to these ideal temperatures by a transfer function G^(s) as follows: 

Tj ( s ) =* G j (s)Y j (s) ( 1 3 A ) 

where s is the Laplace transform variable and j represents any of the 
indices t, s, or h. Furthermore, the transfer functions are assumed to be 
of the form: 


Gj ( s ) 


“St 


O, j 


‘c.J 5 * 1 


(14ft) 


This represents a simple first-order lag with time constant t^ y and dead 
time, t Q . Knowing the actual temperatures, the ideal temperatures can then 
be found' from the expression: 


T,(s) 

Y . ( s ) * (15A) 

3 G (s) 


Since the temperatures are not actually known as continuous functions of 
time, but rather as discrete data points, it is convenient to use a different 
transformation, much like the LaPlace transform, but known as the Z-transform 
and defined: 



-k 


k=0 
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where: 


N number of data points 

T time interval between points 

k an integer 

The points are assumed to be taken at uniform intervals, To determine the 

ideal temperatures, an analogy to equation (15A) is used in the Z-domain: 

Tj ( z ) 

Y. { z ) = (17A) 

3 G (z) 


Transforming equation (16A) into the Z-domain, and making use of the defini- 
tion of the transform in equation (16A) yields the end result: 


W 


T 3 ( fc i) 




1 - e 




(18A) 


where t ± represents the time at which the i th temperature reading was taken. 
It can be shown equation (18A) is valid even if t - t i is not equal to 

T. 


The various time constants, t . and dead time constants, t . are 

C f j O f J 

determined experimentally by running the rig until the temperatures have 
reached a steady state, then creating a sudden change in the load applied to 
the bearing. A time constant is simply the time it takes for a temperature to 
reach 63*2 percent of its new steady state value, while the dead time constant 
is the time required for the temperature to begin responding after the appli- 
cation of the load step. In this case, t , and t are found to be 
approximately 550 sec and t^ to be 700 'sec, while the dead time con- 
stants are substantially less' than T, and thus negligible. 

Equation (18A) can be applied to the actual temperature data to obtain 
the temperatures that would have resulted if the system had negligible thermal 
inertia. These ideal temperatures can be substituted into equation (10A) in 
place of the actual temperatures to obtain a test bearing torque measurement 
with a much more rapid response time. 

There is one notable drawback to this method of lag removal. Any noise 

/( ^cj) 

in the data gets multiplied by a factor of 1/(1 - e J. For a nominal 

sampling time of T = 20 sec and a time constant of t c = 600 sec, this 
means approximately a 30-fold increase in the noise. The torque values 
calculated from the resulting idealized temperature data tend to be so noisy 
that small effects are obliterated. This method of lag removal still has some 
merit for larger effects, however, and the experimentalist is left to choose 
between slow response and substantial noise in this effort to interpret the 
data. 
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TABLE I. - PHYSICAL PROPERTIES OF TEST LUBRICANT (SRG 200) 

[Reference 14.] 

Kinematic viscosity, CS 

at 38 ’C 750 

at 99 °C 42.2 

Viscosity index . . . . . ... 108 

Pour point, °C -7 

Flash point, °C 288 


TABLE II. - BEARING GEOMETRY AND OPERATING PARAMETERS 

Bore, m 0.04 

Ball diameter, m 7.9375xl0" 3 

Pitch diameter, m 5.40xl0~ 2 

Contact angle, deg 12.7 

Ball complement 21 

Nominal load, N . . ........ 891 

Mean outer race stress, N/m 2 . l.OxlO 9 

Mean inner race stress, N/m 2 1.2xl0 9 

Contact dimensions 

Outer semimajor axis, m 4.98xl0” 4 

Outer semiminor axis, m ....... 1.22x1 0“ 4 

Inner semimajor axis, m 5.13xl0” 4 

Inner semiminor axis, m 1.04xl0' 4 
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Figure 1. — Overall experimental apparatus. 
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Figure 2.—' TEHL apparatus. 
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Figure 5. — Electrical connections-LCZ meter to test bearing. 
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Figure 7. — Heat flow in outer housing. 


19 



Capacitance, pF Bearing conductance, ps « micro Siemens 


2.6 




0 5000 10 000 15 000 


Elapsed time, sec 

(b) Capacitance across the bearing versus elapsed time. 
Figure 8.— Electrical measurements. 
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(b) Torque output of outer race drive motor versus elapsed time. 
Figure 10.— Temperature and torque measurements. 
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(c) Bearing torques calculated from motor torque and temperature data. 
Figure 10. — Concluded. 
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Figure A1. — Electrical analogy to heat transfer in bearing housing. 
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